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An integrated marine atmospheric boundary layer (MABL) 
model is evaluated relative to its ability to estimate long- 
wave and shortwave radiation. The model is initialized and 
verified using data taken during the 1983 Mixed Layer 
Dynamics Experiment (MILDEX). Model computations of short- 
wave and longwave radiation are compared with measurements 
made during both atmospheric frontal and non-frontal situ- 
ations. The model results did not always agree with -the 
measurements but reasons for them seem to be known and are 
discussed. One problem is that the MABL model only predicts 
clouds in the boundary layer and does not consider upper 
clouds. This led to most of the major differences. Further 
development of the model with upper-layer cloud specifica- 
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fe oN Repvueril© 


Understanding properties of and processes which affect 
the atmospheric boundary layer is crucial to Department of 
Defense operations. For example, the performance of nearly 
all electro-magnetic (EM) and electro-optical (EO) systems 
is affected by conditions in the marine atmospheric boundary 
layer (MABL). Signal performance is distorted by refrac- 
tion, wave front distortion and extinction. Refraction, 
which leads to the formation of ducts and the trapping of EM 
Signals, is determined by the vertical gradient of the index 
of refraction, which is a function of the vertical gradients 
of temperature, humidity and pressure. EO systems are 
affected by small scale inhomogeneities in the index of 
refraction (due to turbulence), and by water vapor and 
aerosol concentrations in the MABL. Extinction is due to 
water vapor absorption and marine aerosol scattering in the 
presence of-high humidity, and thus is normally restricted 
to the MABL. 

Radiation fluxes coupled with atmospheric dynamics and 
thermodynamics is recognized as one of the major features of 
the MABL. Clouds are the most dominant components in 
affecting the radiation budget of the atmosphere and in 
modifying heating and cooling features of the surface and at 
eoud tops. 

The height of the MABL is controlled to a large extent 
by thermal forcing. During the day, heat is added to the 
surface and is transported upward by buoyant thermal turbu- 
lent eddies. The impingement of the stable layer by these 
eddies entrains overlying air so that in the absence of 
subsidence, the inversion height will rise as long as there 
is an upward heat flux at the surface. The turbulence of 


the cloud topped PBL is also driven by radiational cooling 


at the cloud) top: Hence, radiational cooling at the top of 
the mixed layer is effective at producing entrainment. 

The purpose of this thesis is €o6 examine Ene abil ee 
estimate radiation in the boundary layer with a simple MABL 
model. The model used is a MABL model which has been devel- 
oped at the Naval Postgraduate School (NPS) (Davidson et al, 
1984). This examination was based on the use of simple 
mixed layer parameterizing of specific humidity and poten- 
tial temperature, and the use of a simplified boundary layer 
radiation model. This is in comparison to other schemes 
which could have been used. These include the use of a 
Simplified radiation scheme with detailed atmospheric 
description or the use of a detailed radiation model with 
Simple mixed layer parameters for describing the MABL. The 
scheme chosen may confuse the causes of resulting inaccura- 
cies, since it is not known whether the inaccuracies result 
from incomplete radiation calculations or from incomplete 
vertical profiles. 

The use of potential temperature and specific humidity 
within the well-mixed region and their gradients at the base 
of the inversion also makes this model applicable to stud- 
ying inversion changes. The model diagnoses time evolutions 
of the mixed layer values of inversion height, potential 
temperature and specific humidity, as well as their jumps at 
the inversion (Davidson et al,1984). 

The data used in this study were collected during the 
Office of Naval Research sponsored Mixed Layer Dynamics 
Experiment (MILDEX) between 24 October 1983 and 10 November 
1983. MILDEX is an offshore experiment designed to collect 
intensive Ocean Boundary Layer (OBL) and MABL data, for 
model verification. MILDEX was conducted in the eastern 
north Pacific Ocean, 150 miles off the coast of Lompoc, 


California. 
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A secondary objective of this thesis is to present a 
thorough description of synoptic features occurring during 
MILDEX. Therefore, synoptic descriptions beyond that needed 
for radiation-model considerations are presented in these 


discussions. 


Hal 


Ii. MODEL DESCRIP ECN 


Boundary layer research has led to several time depen- 
dent MABL models based on entrainment energetics and radia- 
tive fluxes (e.g. Deardorff, 1976; Stage and Businger, 1981; 
Davidson et al., 1984). 

The Naval Postgraduate School (NPS) model (Davidson et 
al, 1984), used in this thesis, is a zero-order two-layer, 
integrated mixed-layer model, consisting of a well-mixed 
turbulent boundary layer underneath a relatively non- 
turbulent free atmosphere separated by an inversion (or 
"transition zone"). In a zero-order model, the inversion is 
assumed to have zero thickness and hence profiles of conser- 
vative variables show a discontinuity or a "jump" at the 
inversion rather than a finite gradient (see Fig. 2.1). 

In this thesis, the model was used as a diagnostic to 
obtain radiation, given vertical profiles of potential 
temperature and specific humidity. It is noted that the 
model also can be used to predict profiles and cloud condi- 
tions, and hence radiation at future times. 

Total specific humidity, q,, and equivalent potential 


temperature, O are conserved quantities in  pseudo- 


e/ 
adiabatic processes and are thus assumed to be well-mixed 
within the boundary layer. They are given by the following 


equations: 

Ge oda no (2.1) 
where Slog is the water vapor mixing ratio and qy is the 
liquid water mixing ratio computed using procedures by 


Deardorff (1976); 


Op es ey Ch) (2225 
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where © is the potential temperature, L,, is the latent heat 


Vv 
of vaporization of water, Cp is the specific heat at 
constant pressure. A third fundamental dependent variable 


1s the mean mixed layer depth, h. 

Turbulence is responsible for the transport of kinetic 
energy, sensible heat and latent heat at the level of the 
inversion. Hence, turbulence is responsible drag the 
coupling between the free atmosphere and the MAEL. The 
turbulent fluxes of the temperature and humidity at the 
surface and the level of the inversion change the values of 
the well-mixed quantities over time. Large velocity fluctu- 
ations and mixing in the ABL lead to energetic eddies 
extending from the surface to the inversion. These eddies 
entrain warm, dry air into the mixed layer resulting in 
upward growth of the mixed layer. iis ~ 1s * <called 


entrainment. 


ue 


Another important property of the MABL is the lifting 
condensation level (LCL), which is a function of temperature 
and humidity values in the mixed layer. The LCL determines 
the height at which moisture condenses within an air parcel 
that is lifted adiabatically. The height of the LCL rela- 
tive to the height of the inversion specifies cloud forma- 
tion and dissipation within the mixed layer. If the LCD wie 
above the level of the inversion, the mixed layer is cloud 
free. If the LCL is below the inversion, a layer of stratus 
is formed, extending from the LCL to the top of the mixed 
layer. Stratus in the mixed layer has profound effects on 


the radiation budget. 


A. RADIATIVE FLUX CONSIDERATIONS 

Radiation, as a basic driving force of the diabatic PBL, 
is an important component of the surface energy budget. 
Cloud-radiation effects have a significant influence on 
vertical cloud development since radiational cloud top 
cooling increases convective activity (Beniston and Schmetz, 
IIe >). ; 

The following section describes radiation scaling in the 
model. It 1s an abbreviated version of that given by 
Davidson et al (1984). The cloud-shortwave radiation inter- 
actions are modeled on the assumption of solid stratocumulus 
© loud seover. The longwave model is based on an empirical 
relation of the average cloud liquid water and the cloud 
emissivity which may contain some slightly broken thin stra- 
tocumulus. Clearly neither radiation model applies under a 
field of individual optically thick clogs. 

1. Lonqwave Radiative Flux 

The longwave net radiation is the amount of 
upwelling longwave radiation measured at the surface which 
includes a positive contribution from surface emittance up 
and negative (downward) contributions from the cloud base 


and molecules including water vapor. Longwave radiation is 
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a function of sky temperature, water content, and cloud 
emissivity. The long-wave radiation flux was modified from 
a scheme used by Stage (1979), to permit non-black stratus 


clouds by introducing the emissivity, € which is a func- 


Gc 
tion of the total liquid water content, W, 


i. = Ll = exp (aW) =e C23) 


and 


(2.4) 
W= | (p q,) da, 
Zo 


where p is the density of air, q, is the cloud liquid water 
specific humidity, Z is the height, ae is the LCL (cloud 
bottom) and a = 0.158 m°g * (Slingo et al, 1982). Since the 
cloud liquid water profiles are approximately linear with 


height, Eq.2.4 becomes 
ie 8-5 9 (h = 2) Gih » (GZ >) 


where Qj, is the liquid water content at the cloud top. 
The long-wave cloud top net radiation flux, Bee la: 
calculated from the cloud top temperature, oe and the 


effective radiative sky temperature, ak using the 


sky, 
Stefan-Boltzman law 


path SRS Ca ae ae (2.6) 


where ¢ is Stefan's constant and &. depends on the inte- 


grated liquid water content (see Eq. 2.3). Similarly, the 


mom at the cloud bottom, er is 


= 4 4 
ie = cao (T.- = To), C257) 


ILS 


where T, is the sea-surface temperature and T. is the cloud 
base temperature. For the cloudy case, the model neglects 
the flux divergence in the clear air between the sea surface 
and the cloud bottom, therefore, Li, = Ly, 
net longwave radiation at the surface. 


where Ee is the 


Distinctions must be made for cloud-free and cloudy 
cases. In the cloud-free case, the net longwave flux is 
calculated from the water vapor and temperature profiles. 
The net fluxes By and ae are calculated at Z=h and Z=0, 
respectively, using procedures described by Fleagle and 
Busrnger "(1930 )- These procedures include integrating the 
flux emissivity profile with the previously mentioned modi- 
El1Gacaon, JEQs w2so. In the cloudy case, only the 24> h 
part of the integration is required in order to obtain the 
effective sky temperature. 

The simplified longwave radiation equations treat 
the boundary layer as a two-layer model. The effects of 
continuum absorption are not addressed by the _ model. 
Continuum absorption by both water dimer and water vapor is 
important in the lower part (4 km) of the atmosphere ie 
moisture values are large (Selby et al, 1976). This leads 
to a Known bias in the longwave radiation calculations. 
This bias is to reduce downward emitted longwave radiation. 

Zz. Shortwave Radiative Flux 

The model calculates the shortwave radiation from 
the solar angle, which is a function of latitude, julian day 
and the time of day, and meteorological parameters. The 
meteorological parameters include liquid water mixing ratio, 
total mixing ratio, inversion height and mixed layer wind 
speed. The wind speed is used to estimate local production 
of aerosols. In calculating downward solar shortwave 
radiation received in the surface layer, the model 
calculates downward direct shortwave radiation and the 
downward diffuse shortwave radiation scattered by cloud 


droplets and aerosols. 
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The following shortwave radiative flux discussion is 
brief. The complete discussion and equations used in the 
model are given by Davidson et al (1984). A delta Eddington 
shortwave radiation flux calculation (Joseph et al,1976) is 
used in the model. The solar flux is evaluated over 15 
equally spaced (0.1 wm width) bands from 0.2 [lim to 1.7 Im 
wavelength. The incidence flux at the top of the mixed 
layer is obtained from the flux at the top of the atmosphere 
and the average transmitance in each of the fifteen bands, 
using the data and methods in the Smithsonian Meteorological 
Tables (1963). The fact that the solar flux in the model is 
only evaluated from absorption bands less than 1.7 [Itm wavel- 
ength, is a deficiency since a major water vapor absorption 
band is at 1. 87 [|tm. 

Absorption is due to five water vapor absorption 
bands centered at 0.75, 0.94, 1.10, 1.38, and 1.87 itm wavel- 
engths. Using the water vapor absorption calculated from 
LOWIRAN 3B (Selby et al, 1976), the average absorption coef- 
ficients in each of the relevant O.1] tm width bands (0.7 to 
1.7) were calculated as a function of total absorber amount 
(precipitable water) (Davidson et al,1984). 

In addition to water vapor absorption, there is 
scattering and absorption in the mixed layer by atmospheric 
particles, cloud droplets in the cloudy case and sea salt 
aerosols in the clear case. The short wave calculation is 
restricted to the region Ze Jah In the cloud-free 
case, ae is set to zero to include the short-wave attenua- 
tion in the entire (cloud-free) mixed layer. In the cloudy 
case, only the region from the cloud base up to the inver- 
sion is considered in the calculation, because aerosol scat- 
tering below the cloud base is much less than scattering 
above the cloud base. 

The scattering properties of the particle are calcu- 


lated using the Mie coefficient approximations, given by 


iy) 


Deirmendjian (1969), at each band wavelength for the speci- 
fied particle size distribution in each layer. Partieme 
absorption is accounted for by using the complex refractive 
index for pure liquid water from Hale and Querry (1973). 
The absorption effect of water vapor and scattering plus 
absorption effects of the particles are combined in the 
manner suggested by Lacis and Hansen (1974). The particle 
size spectra are specified as log-normal distributions and 
the scattering calculations are done for four size intervals 
equally spaced in log-radius space. 

The sea-salt-aerosol parametization is an empirical 
fit to data published by Fairall et al (1983). The cloud 
droplet spectrum is based on empirical parametizations of 
field measurements. The data indicate a relatively constant 
total number density and a steady increase in mode radius 
With height above the cloud base. This 1s because the 
liquid water content, which describes the droplet size 
distribution, depends on the height above the cloud base 
(LCL). Hence, shortwave cloud transmittance in the model is 
dependent on liquid water content (Slingo and Schrecker, 
TIo2) 3 
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TIL. MILDEX DATA COLLECTION AND SYNOPTIC DISCcUssloO 


A. MIXED LAYER DYNAMICS EXPERIMENT (MILDEX) 

The Mixed Layer Dynamics Experiment (MILDEX) is a multi- 
group and multi-platform experiment. The observations were 
taken aboard the R/V Acania as well two other platforms, the 
R/P Flip and R/V Wacomo. The observations were made off the 
central coast of California near 34° N, 124° W between 24 
October 1983 and 10 November 1983 (see Fig. 3.1). 

The purposes of MILDEX were: 


‘LO Be ovads magnitudes of air-sea exchange rates for use 
as boundary conditions in mixed layer modeling; 


¢ to evaluate the drag coefficient from turbulent kinetic 
ore dissipation measurements and its dependence on 
swell amplitudes and direction; 


¢ to provide a time series of MABL structure from radio- 
sonde measurements for model verification; 


¢* to provide measured radiative data for comparison with 
radiative transfe r algorithms developed at Scripps 
Institution of Oceanography; 


¢ to evaluate the effectiveness of SODAR as a ship-borne 
instrument to measure changes in MABL parameters; 


¢ to evaluate radiative transfer models of the atmosphere 
with cloudiness as a primary parameter. 


Measurements made from the R/V Acania included those of 
windspeed, temperature, humidity, radiation and atmospheric 
pressure in the surface layer, and sea-surface temperature. 
Profiles of wind speed and direction, temperatures and 
humidity were obtained from on-board radiosonde launches. 
Radiosonde launches were coordinated with satellite pass 
times. Additionally, there were hourly observations of 
cloudiness and sea _ state. Figure 3.2 is a time series 
depicting the wind speed, temperature and humidity for the 
entire MILDEX period. A time series of net shortwave and 


net longwave radiation at the surface is shown in Fig. 3.3. 
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Peover Llic CONDITIONS 

This section describes atmospheric conditions during 
MILDEX. Emphasis is placed on conditions appropriate to the 
mixed layer model and radiation variations. These condi- 
tions include position of fronts, cloud conditions and what 
clouds were observed in and above the mixed layer. Five 
frontal passages occurred during the MILDEX experiment. 
Major frontal passages occurred on 29 and 31 October, and 6 
November 1983. 

In the period from 25 October to 27 October, synoptic 
patterns were dominated by a subtropical high pressure cell. 
Winds were light (1-6 m/s) and northerly. Dominant clouds 
were stratus, stratocumulus and cirrus. The top of the 
inversion was near 700 m, and potential temperature and 
specific humidity below the inversion remained well-mixed 
throughout this period, 

On October 28, a weak cold front began to influence 
conditions. By l1800GMT winds shifted to southeasterly and 
wind speeds decreased from 5-6 m/s to 1-3 m/s. Air tempera- 
ure increased due to an advancing warm air mass. Croud 
cover was mostly broken to overcast dominated by stratocu- 
mulus and altostratus (see Figs. 3.9, 3.10 and 3.12). 

The front passed at I1400GMT on 29 October, even though 
it was not shown on the National Meteorological Center (NMC) 
surface analysis in Fig 3.14b. The 500 mb analysis showed 
Mmimes of a 500 mb trough near 130° W, with a low at the top 
of the trough at 30° N. The frontal passage was determined 
by satellite interpretation and a time series analysis of 
observed winds (speed and direction), observed temperatures 
and pressures (see Fig. 3.2). Winds shifted from southeast- 
erly to southerly, and wind speeds increased from 3 m/s to 
12 m/s before the frontal passage and then decreased to 2 
m/s after frontal passage. The warm sector preceding a 


second cold front increased the temperature by two degrees. 
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Sebo eames es: 26 October 1983 Surface Charts 
a) 1200GMT (top), b) 1800GMT (bottom). 
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a) OOOOGMT (top), b) O600cMT (bottom) 
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Fig. 3.10 28/2045GMT October 1983 Satellite Picture. 
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sec. Sashes 28 October 1983 Surface Charts 
a) OOOOGMT (top), b) O600GMT (bottom). 
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Cloud cover ranged from scattered to overcast. Cumulonimbus 
clouds were reported at 1500GMT, along with rain from 
O6OO0GMT to 1400GMT (see Fig. 3.13). 

On 30 October 1983, the next cold front remained west of 
the experiment region. Sum Pace pressures increased, and the 
wind shifted from southerly to westerly with wind speed 
steadily decreasing from 8 m/s to 2 m/s. Air temperatures 
fell one degree. The cloud cover was clear to scattered 
with cumulus and stratocumulus clouds present (see Fig. 
Seloa ). 

A weak cold front passage was reported by the R/V Acania 
at about 1700GMT on 31 October (see Figs. 3.15b, and 3.16). 
The synoptic picture shows a low in the eastern north 
Pacific Ocean, and a high dominating the Pacific north of 
Hawaii. The entire period was characterized by a 100% 
ceiling and intermittent rain and drizzle. Several squall 
lines were noted. after the frontal passage. The clouds were 
mainly altostratus before the frontal passage and stratus 
after the frontal passage. . The winds shifted from south 
southwest to west after frontal passage. The wind speed was 
7 to 9 m/s ahead of the front and 3 to 5 m/s behind it. The 
air temperature dropped by 1.5 degrees centigrade after the 
front went by. 

Tieceeoassage Of thes third cold front occurred on 1 
November 1983 at approximately 1200GMT, accompanied by a 
deep 500 mb trough with a northwest to southeast tilt (see 
Peeas. 3.17 and 3.18). The winds shifted from westerly to 
northwesterly with wind speeds increasing to 12 m/s at 
16Q0O0GMT. During a 30-minute period, the temperature dropped 
by 1.5 degrees centigrade. Over the next two days, the 
subtropical high pressure cell re-established itself. Wind 
directions were northerly to northwesterly, at speeds of 2 
to 5 m/s. The temperature began to drop on November 3, with 


patches of fog, rain, and drizzle reported throughout the 
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day, as a stronger cold front was approaching (see Fig. 
Bre) 

The synoptic pattern for 4 November 1983 showed a low 
center in the Gulf of Alaska, and a high building west of 
Baja California and south of the experimental region . No 
frontal passage was observed. Winds were generally from the 
north shifting to the west, then returning to the north, 
while wind speed fluctuated from 3 to 7 m/s. The surface 
pressure and the air temperature were nearly steady. As 
seen on the satellite pictures (see Figs. 3.20, and 3.21), 
it was overcast in the late afternoon on 3 November, with 
total clearing at night, and overcast again by mid-morning 
of the 4th. Fog, cumulus, and stratus were observed in the 
afternoon of the 4th (see Fig. 3.23). 

The experimental region came under the influence of a 
weak pressure cell on 5 November, located to the _ south. 
Another low came in northwest of the experimental site, with 
an associated cold front. By 6 November, the warm sector of 
the fourth cold front began to influence synoptic conditie@> 
(see Figs. 3:26, andes z sa). Winds, which were southwes-= 
terly, increased in speed from 2 to 13 m/s. The advancing 
warm air mass increased air temperature by 1.2 degrees 
centigrade. Cloud cover was broken most of the day, with 
cumulus and stratus increasing during the evening due to the 
approaching front. Sea level pressure increased to a 
maximum after the cold front passed at OOOOGMT on the 7th of 
November (see Fig. 3.28b). The wind shifted to northerly 
with wind speeds averaging 10 m/s. A 500 mb ridge followed 
the front and was over the experimental site at approxi- 
mately 1Z2OOGMT on the 8th (see Figs. 3.29, and 3.3la). 
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a) 1200GMT (top), b) 1800GMT (bottom). 
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IV. RESUS 


A. MODEL APPLICATION 

Thirty-six MILDEX radiosonde soundings were examined for 
possible radiation diagnosis using the MABL model. Those 
examined were based on the potential temperature and the 
specific humidity profiles (see sample output Fig. 4.1). 
All cases with a distinct humidity inversion were further 
examined relative to radiation. The inversion height, the 
mixed layer values oF potential temperature and specific 
humidity, the potential temperature and specific humidity 
jumps across the inversion, the potential temperature and 
specific humidity lapse rates above the inversion, and the 
LCL were calculated from the profiles. These calculated 
values were used for a diagnosis of radiation. Also used 
were the Julian day, time of radiosonde launch, latitude, 
surface pressure and sea-surface temperature. 

In this study, the interest was on how an integrated 
model using a simplified radiation scheme diagnosed radia- 
tion at a specific time. Therefore, only results from the 
first time step in the model calculations were used in the 


evaluation. 


B. EXPLANATION OF TABLE I 

Table I contains model and observed results for the 
selected soundings. FL refers to the radiosonde launch 
number. Of 36 launches, only the 17 listed had the features 
meeting the mixed layer model (see Fig 2.1). The date is in 
Pacific Standard Time (PST). Net incoming shortwave and net 
outgoing longwave radiations measured at the surface are in 
watts per meter squared. Measurement errors (uncertainties) 
for shortwave values are estimated to be 3% or nominally +2 
W/m? from 15 to 60 W/m* and +10 W/m* for midday values around 
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Fig. 4.1 Sample Model Potential Temperature and 
Spectitecehumztailty Profiles = valid 27/03452. 
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300 W/m*. Measurement error and surface estimated error for 
longwave values are +15 W/m? (personal conversation, K. 
Katsaros, Atmospheric Science Department, University of 
Washington at Seattle). Clouds types are abbreviated and 
are in tens of % (1.e. 2 cu is 20% cumulus). N/7Geane@e2eeue 
29 shortwave radiation column refers to not considered due 
to measurement error. 

The model clouds in table I are based on model diagnosis 
of the presence or absence of a cloud within the PBL. The 
model uses the LCL as the cloud base, and computes the cloud 
depth by subtracting: the LCL from the height of the inver- 
sion. If the LCL is higher than the inversion, there are no 
clouds in the boundary layer, and the cloud depth is set to 
zero. The observed clouds in table I were based on both 


in-situ observations and satellite imagery. 


Gy “DISCUSSION CE RESULTS 

in? “this section, the similarities and differences 
between the model and the observed results are discussed. 
On the morning of October 26, a high pressure system domi- 
nated, and the observed sky was 20% cumulus (see Figs. 3.4 
and -335:)2 The model LCL was 157 m below the inversion, 
leading to a model cloud that was 157 m thick. The model 


cloud cover was, of course, overcast and extended from 


horizon. “Com neraezon. The model shortwave radiation (60.2 
W/m?) from 10/26/0800 was 20% greater than the hourly aver- 
aged measurement (50.3 W/m°*). This discrepancy occurred 


around sunrise when the solar angle was low, and accurate 
measurements are diftrrveult and time sensitive. 
Additionally, the observed values were an hourly average 
which includes values taken before and after O800 PST. The 
model outgoing (net) longwave radiation (36.6 W/m*) was less 
than the observed value (71.8 w/m*). The model diagnosed a 
relatively thick (157 m) overcace Clouds emit longwave 


radiation back to earth. Hence the absence of observed 
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clouds caused the larger observed net longwave radiation to 
space. 

Later that night (10/26/2100), the high still dominated 
(see Fig. 3.7) but the sky became overcast ( 100% SeDR The 
model diagnosed a thick cloud (170 mm), and since the 
observed cloud condition was overcast, the longwave radia- 
tion values were nearly equal (model: 37.2 , obsvd: 37.8 
W/m°). Since 2100 PST was well after sunset, short wave 
radiation was not considered. 

The high remained the dominant synoptic feature for the 
next period (10/27/1/00) (see Fig: “sai The observed 
cloud condition was overcast (100% sc) but the satellite 
picture (see Fig. 3.8) showed breaks in the overcast. The 
short wave radiation value diagnosed by the model (56.9 
W/m?) was almost equal to the measured value (55.4 W/m’). 
This 1S a surprising result since the model diagnosed no 
clouds, and the observed stratocumulus deck was completely 
overcast. This period is just before sunset, so the zenith 
angle of the sun was low. It is unknown if the horizon was 
actually clear, which would have explained the model's good 
performance. The longwave diagnosed value (191.2 W/m*) was 
over six times the measured value (29.7 W/m*). The differ- 
ence is due to the fact that the observed clouds emitted 
longwave radiation back to earth and decreased the observed 
net longwave radiation. 

Through the next period (10/28/0630) stratus cloud 
coverage existed (see Fig. 3.9) and the high pressure system 


was still still the dominant synoptic feature (see Fig. 


ee WA Since this time was before sunrise, shortwave radia- 
tion was not considered. The model diagnosed a thin cloud 
(74 m), which emitted longwave radiation back to the 


surface, but not as much as the observed full dense stratus 
deck. Therefore the model overestimated the net longwave 
radiation (63.9 vs 23.0 W/m‘). 
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Later on the 28th (1730), the model still underestimated 
cloud coverage. A broken (60%) lower layer of cumulus was 
Beservea dn ada@it10on “to high clouds (see Fig: 3.15). The 
high began to weaken as the warm sector ahead of a cold 
front approached (see Fig. 3.14a). The model diagnosed LCL 
was close to the inversion (38 m below), so a thin cloud was 
estimated. The high clouds emitted more longwave radiation 
back to the surface than did the model's environment of no 
clouds above the boundary layer (model: 108.2, obsvd: 25.8 
W/m*). The model's shortwave diagnosis (32.3 W/m*) was 
consistent with the- observed value (28.0 W/m*) since the 
model diagnosed a thin cloud and broken clouds were 
observed. The concern with this case is the agreement in 
the shortwave radiation values but large disagreement in the 
longwave values results. Perhaps, the shortwave radiation 
agreement exists because of a low sun angle. 

For =the next time period (10/2971700), the cold front 
remained west of the experimental region (see Fig. 3.15a). 
Observed cloud cover was scattered cumulus and stratocumulus 
clouds. Although the model diagnosed no clouds, the model 
still underestimated the measured shortwave radiation 
(model: 55.4; obsvd: 85.6 W/m*). This period is near sunset, 
a time not optimal for testing of the model. Observed upper 
clouds in addition to scattered cumulus emitted longwave 
radiation which, caused the model (182.9 W/m*) to overesti- 
mate the net longwave radiation observed (55.9 W/m’). 

The next period (11/1/0530) coincided with a frontal 
passage (see Fig. 3.18a). Cloud cover increased (see Fig. 
3.17) and was accompanied with rain, showers and drizzle. 
The valid time was before sunrise, thus shortwave radiation 
was not considered. The model again failed to predict any 
clouds, however clouds were observed above the boundary 
layer (50% cc). Therefore, the model, without upper clouds, 
overestimated the net longwave radiative flux (197.6 vs 44.9 
W/m). 
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The subtropical high was re-established on 11/3/1730 
over the MILDEX area (see Fig. 3.23a). An overcast stratus 
deck was observed (see Fig. SNZOn x The time was near 
sunset (1730) and the model's shortwave diagnosis (23.4 
W/m*) agreed with the observed value (26.9 W/m*) since the 
model diagnosed a cloud and stratus clouds were observed. 
The model LCL was 39 m below the inversion, producing a 
relatively thin cloud. The observed stratus deck emitted 
more longwave radiation than the 39 m diagnosed cloud, so 
the model overestimated the net longwave radiation (141.9 vs 
9.2 W/m). 

A high pressure system continued to dominate the next 
period (11/5/0930) examined (see Fig. 3.2565). Scattered 
stratus and broken cirrus were observed and the model diag- 
nosed the occurrence of thin boundary layer clouds (28 m) 
correctly. The model diagnosed the highest value (391 W/m*) 
of shortwave radiation in this data set, and the measured 
value was also one of the highest (305 W/m). The broken 
cirrus clouds reduced the incoming solar radiation, causing 
a higher model diagnosed value than what was actually meas- 
Ured. This case (FL26) was similar to other cases (.e.g. 
FL's 5, 6, 22, and 35), in which the LOB was close tompme 
below the inversion. However, high clouds were observed and 
presumably, these high clouds and the scattered stratus 
clouds emitted more longwave radiation to the surface than 
did the thin diagnosed boundary layer cloud. Thus, the 
model's net longwave emittance was larger (model: 226.7, 
obsvd: 77.5 W/m*). 

A weak cold front to the west dissipated, and had little 
effect on the synoptic conditions during the next period 
(11/5/1730) (see Fig 3.26a) examined. Scattered cirrus 
clouds were the only clouds observed and, as with other 
cases (e.g. FL's 4, 7, 16, 30 and 36), the LCL was diagnosed 


above the inversion for FL 27. However, the observed high 
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clouds emitted longwave radiation back to the surface. 
Thus, the model's clear sky had a larger net longwave radia- 
tion value (model: 184.7, obsvd: 70 W/m*). The model's 
sunset diagnosis of shortwave radiation (35.8 W/m’) was 
again smaller than the measured value (55.7 W/m*) because of 
the diagnosed and observed cloud differences. 

The warm sector ahead of the next approaching cold front 
began to influence synoptic conditions during the next 
periods (11/6/0330 and 11/6/1100) (see Figs. 3.26b and 
3.28a). Cloud cover was observed as scattered to broken 
cumulus throughout the day. FL 28 was a morning launch so 
shortwave radiation was not considered. The model diagnosed 
a low LCL for both periods in the morning of November 6 
fer /6/0330 and 11/671100), resulting in thick clouds (233m 
and 328m, respectively). During both periods, scattered 
cumulus (less than 40%) were observed. In both cases, the 
model underestimated the net outward longwave radiation 
(31.9 vs 46.6, and 35.9 vs 47.6 W/m’, respectively). Since 
the model diagnosed thick clouds extended from horizon to 
horizon, the model results indicated more longwave radiation 
emittance back to the surface than that emitted from the 
observed cumulus. 

The next period (11/6/1730) coincided with a frontal 
passage (see Fig. 3.28b). The observed sky was overcast 
cumulus, but the model diagnosed no clouds. As in previous 
cases (e.g. FL's 4, 7, 16, and 27), the model overestimated 
the net longwave radiation emitted (model:173.3; obsvd: 8.2 
W/m). This period was again near sunset, so the model's 
shortwave radiation diagnosis (28.4 W/m*) was not in agree- 
ment with the measured shortwave radiation (41.8 W/m’). 

During the next period (11/8/0500), an intense low pres- 
sure system was located to the north in the Gulf of Alaska, 
but a weakening high was still the dominant synoptic feature 


in the MILDEX region (see Fig. 3.31la). Stratus clouds were 
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observed (see Fig. 3.29). As with previous cases (e.g. FL's 
1, 28, and 29), the model diagnosed a thick cloud (164 m) 
for FL 33, while the observation was only scattered (20% 
st), therefore the model underestimated net longwave emit- 
tance (71.1 vs 99.8 W/m’). The launch was well before 
sunrise, thus shortwave radiation was not considered. 

Cloud cover increased and became overcast over the next 
period (11/8/1200), (see Fig. 3.30). The high dissipated, 
as the warm front moved into the region (see Fig. 3.31b). 
The model diagnosed LCL was 177 meters below the diagnosed 
inversion. As was the case with FL 2, the model value (75.2 
W/m?) agreed closely with the observed value (67.2 W/m’*), 
slightly overpredicting the longwave emittance. This noon 
launch time produced the highest measured value of shortwave 
radiation (369.0 W/m*) received. The model value (273.5 
W/m*) was lower since the model diagnosed a thick horizon to 
horizon cloud, while the observed cloud cover had breaks in 
the overcast. 

During the next period (11/8/1730), the MILDEX region 
was in the warm sector of an approaching cold front (see 
Fig: 3c. Cloud cover was overcast (100% sc) (see Fig. 
SoZ The model diagnosed a thin cloud (40m). As™ im 
previous cases (e.g. FL' 5, 6, 22, and 26) in which thin 
clouds were diagnosed but broken to overcast clouds were 
observed, the model overestimated the net longwave radiation 
(158.1 vs 35.3 W/m*). This period was again near sunset, so 
the model's shortwave radiation diagnosis (16.3 "W/m 
although close to the measured shortwave radiation (12.4 
W/m*), was not of significance. 

The final period (11/8/2130) considered was prior to the 
final frontal passage of the MILDEX experiment (see Fig. 
333)):. Cloud cover was overcast (see Fig. 3.32). The last 
launch (FL 36) was at night, therefore shortwave radiation 


was not considered. This case is similar to previous cases 
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(e.g. EFL's 4, 7, 16, 27, and 30) when no clouds were diag- 
nosed, but upper clouds were observed (100% cc), therefore 
the model greatly overestimated the net longwave radiation 
emitted (161.8 vs 22.4 W/m’). 


oe) 


V. CONCLUSIONS 


Over long periods of time, energy absorbed (shortwave 
radiation) by the MABL is approximately equal to energy 
emitted (longwave radiation), but predicting the perturba- 
tions in radiative flux is important in the evolution of the 
MABL. The perturbations of boundary layer parameters such 
as temperature and humidity can change in a matter of hours, 


thus affecting and being affected by radiation flux. 


A. SHORTWAVE SUMMARY 

The model calculated shortwave radiation values from 
solar angle and MABL parameters. At night, shortwave radia- 
tion was not considered (e.g. FL's 2, 5, 16, 28, 33, and 
36). Discrepancies occurred near sunrise (.e.g. FL 1) gana 
sunsec. (6.9.5 Bis) 7727 sean Oe The reasons for these 
discrepancies could be threefold. First, did the model use 
the proper solar angle to get its shortwave radiation value 
when the solar zenith angle is most critical? Second, was 
the boundary layer cloud cover at the horizon truly repre- 
sentative of the actual cloud cover? Third, the measured 
values were an hourly average. Near sunset/sunrise values 
were based on measurements both before, which gave a non- 
zero/zero contribution, and after which gave a zero/nonzero 
value. Thus, model runs near sunrise and sunset were not 
good for evaluating the shortwave calculations of the model. 

For FL 26, the model overestimated the observed short- 
wave value due to cirrus clouds above the model's environ- 
ment. The model (e.g. FL 34) underestimated the shortwave 
radiation when a thick horizon to horizon cloud was diag- 
nosed by the model and the observed sky conditions indicated 
breaks in the overcast. However, over one half of the 
model's daytime estimations were within 20 W/m* of the meas- 
ured shortwave value (e.g. FL's 4, 6, 22, 27, 30, and 35). 
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Another reason for errors was the simple radiation 
scheme used by the model. Water vapor absorption beyond 1.7 
um was ignored despite the fact that a major absorption band 
exists at 1.87 Um. Also, the shortwave calculation with 
clouds was restricted to the region between the LCL and the 
inversion. In the cloudy case, the layer below the cloud 
base was assumed to emit zero shortwave radiation flux, 


despite the presence of aerosols. 


B. LONGWAVE SUMMARY 

When thick clouds were diagnosed and broken or overcast 
conditions were observed, the model was within 15% of the 
measured value (e.g. FL's 2, and 34). When either no clouds 
or thin clouds were diagnosed, and clouds were observed, the 
model greatly overestimated the net longwave radiation 
Paeeeecmie. Go. FL s 4, 5, 6, 7, 16, 22, 26, 27, 30, 35, and 
36). The model is deficient in that it only diagnoses 
clouds in the boundary layer. 

When a thick cloud (cases with a model cloud greater 
than 157 m) was diagnosed, and the observed sky was scat- 
tered (20% or less), the model underestimated the longwave 
radiation emitted (FL's 1, 28, 29 and 33). The diagnosed 
thick clouds emitted more longwave radiation back to the 


surface than the observed scattered clouds emitted downward. 


C. RECOMMENDATIONS 
These recommendations pertain to specifications within 


the model and to procedures for evaluating it with in-situ 


data. Since the model provides a point estimate, clouds are 
either there or not there (100% or 0%). When clouds are 
diagnosed they extend from horizon to horizon. Partial 


Cloud cover certainly has to be considered for both short- 
wave and longwave radiation effects on the Boundary Layer. 
A statistical approach for percentage cloud cover could, 


perhaps, be included for more realistic cloud effects. 


Oy, 


Improving parametization is necessary for determining 
the LCE: This is important because the height of the LCL 
determines the cloud thickness. Further development of 
mixed layer gradients within the model is necessary. 

The model does not consider clouds above the mixed layer 
and that is a severe limitation. The PBL is not self- 
deterministic, but rather dependent on larger scale external 
factors, such as the radiative effects of upper level 
clouds. 

Many MILDEX radiosondes were launched to coincide with 
satellite passes, and as a result, radiosondes were launched 
near OO000Z and 12002, which corresponds to 0700 PST and 1900 
Boa. Solar shortwave radiation values should be evaluated 
at local noon when the shortwave values are at their 
greatest and the solar angle is least critical. Therefore, 
it is recommended that future radiation data-gathering 
experiments include radiosonde launches close to local noon. 

The model is a potentially powerful tool but still has 
limitations. Curpently, the model is best used for 
explaining dynamical processes of the boundary layer rather 
than radiative processes. Further research and development 
of the model is needed with the eventual goal of an opera- 


tional boundary layer model. 
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